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It is demonstrated, that chirality violating condensates in massless QCD arise en- 
tirely from zero mode solutions of Dirac equations in arbitrary gluon fields. The model 
' is suggested, where the zero mode solutions are the ones for quarks, moving in the 

instanton field. Basing on this model were calculated the quark condensate magnetic 
' c*j ' susceptibilities of dimensions 3(x) and 5 (k and £). The good considence of the values 

Q-f X) K an d £, obtained in this approach with ones, found from the hadronic spectrum 

O ia a serious argument in favour, that instantons are the only source of chirality vio- 

lating condensates in QCD. The temperature dependence of the quark condensate is 
discussed. It is shown that the phase transition, corresponding to the T-dependence 
of the quark condensate a(T) as an order parameter, is of the type of crossover. 
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It is well known, that because of small values of light quark masses the perturbative QCD 
pocesses the property of chiral symmetry in cases, when the heavy quarks contributions can 
^ be neglected. However, in the real hadronic world the chiral symmetry is badly violated. 

This statement evidently follows from the existence of large proton mass and the absence of 
the negative parity baryon, degenerate with proton. It is also well known that large values of 
light quark condensates indicate the violation of chiral symmetry in QCD. These two facts 
are deeply interconnected: the values of the proton mass can be expressed through the value 
of quark condensate [1]. 

In this talk I discuss the appearance of quark condensate in nonperturbative QCD. It 
is established the connection of quark condensate value with zero mode solutions of quark 
Dirac equation [2]. In this way it is clarified the nature of proton mass and, as a consequence, 
the nature of all visible mass in the Universe [3]. 
Consider QCD action in Euclidean space-time 



S= 4 



- J d^xG\ v -J d 4 xJ2 i>f (^V M + im f )ip f 
where is gluon field tensor, the sum is over quark flavours. 



(1) 
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V, = d, + ig^-A; (2) 

and A™ is the gluon field. Pay attention, that in Euclidean formulation of QCD ip is replaced 
by ip + . (The review of Euclidean formulation of QCD and instantons is given in [4], see 
especially [5].) The Dirac equation for massless quark in Euclidean space time has the form: 

-il^^n(x) = Xni'n(x) (3) 

where ip n (x) and X n are the eigenfunctions and eigenvalues of the Dirac operator —V = 
— ry M V M . Expand the quark fields operators into the left and right ones 

^ + = ^4 (1 + 75)+ ^4 (1 " 75) ' (4) 

where 

75^L = IpL, Ib^R = -lj>R (5) 

Then for nonzero A n the Lagrangian and the action reduces to the sum of two terms 



L = 



d 4 x (6) 



completely symmetric under interchange L < — > R. Therefore the solutions of the equations 
for left and right quark fields are also the same - the states, constracted from left and right 
quarks are completely symmetrical. This conclusion was obtained for fixed gluon field. It is 
evident, that the averaging over the gluon fields does not change it. Quite different situation 
arises in case A = 0. The contribution of this term to the Lagrangian: 

AL = J d 4 x[^ + ^]Wo (7) 

is equal to zero and no conclusion can be done about the symmetry of states build from 
left and right quark fields. One of the consequences from the said above is that all chirality 
violating vacuum condensates in QCD arise from zero mode solutions of Dirac equations (3). 
These general arguments are supported by the well known facts: 

1. The general representation of the trace of quark propagator S(x) is expressed through 
the spectral function p(A) as a function of eigenvalues A (Kallen-Lehmann representation): 

TrS(x 2 ) = - [ d\p(\)A(x 2 ,\) (8) 

7T J 

At x 2 = A(x 2 , A) reduces to 5(A) and we have (in Minkowski space-time): 

p(0) = -7T (0 | ^(0)^(0) | 0). (9) 
(The Banks-Casher relation [6]). 

2. The zero-mode solution of (3) for massless quark in the instanton field is the right 
wave function - i/jr(x) = (1 — 75)^(2;) and in the field of anti-instanton is the left one - 

i> L {x) = {l + lb )4>{x) [7,8]. 
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Basing on the statements, presented above, let us formulate the model for calculation 
of chirality violating vacuum condensates in QCD. Suppose, that vacuum expectation value 
(v.e.v.) of the chirality violating operator O c , v , is proportional to matrix element ipoO c , Vi ipo, 
where ip Q is the zero-mode solution of Eq. (3) in Euclidean space-time: 

(0 | tfOcv.ij; | 0) ~ ^ + O c ,>o. (10) 

■00 depends on x, on the position of the center of the solution x c , as well as on its size p: 
■00 — ipo{.x — x c1 p). Eq.(10) must be integrated over x c , what is equivalent to integration over 
x — x c . (In what follows the notation x will be used for x — x c .) We assume, that p =Const 
and find its value from comparison with the known v.e.v. 's. Finally, introduce in (10) the 
coefficient of proportionality n. So, our assumption has the form: 

(0 | ^(0)0,^(0) \0) = -nJ d 4 xij+(x,p)O c . v .Mx,P) (11) 

Our model is similar to delute instanton gas model [9] , where x c is the position of instanton 
center. Unlike the latter, where the instanton density has dimension 4, n has dimension 3 
and may be interpreted as the density of zero-modes centers in 3-dimension space. Note, 
that the left-hand side of (11) is written in the Minkowski space-time, while the right-hand 
side in Euclidean ones. (The sign minus is put in order to have n positive.) In perturbative 
calculation in case, when the number of flavours is more then 1, Nf > 1, the contribution of 
instantons to the action is suppressed by the factor, proportional to the product of Nf — 1 
light quark masses. Therefore, the assumption (11) means that such suppression is absent 
in nonperturbative calculation. For x and p-dependens of ip (x,p) we take the form of the 
zero-mode solution in the field of instanton in SU (2) colour group: 

MX,P) = 2 (l ~ 75 ^ ( X 2 + P p 2)3/2 *0, (12) 

where xo is the spin-colour isospin (| T |= 1/2) wave function, corresponding to the total 
spin I + T = J equal to zero, J = 0. ipo(x, p) is normalized to 1: 

[ d 4 xip + {x,p)ip{x,p) = 1 (13) 

Consider first the quark condensate (0 | qq | 0), the most important chirality violating 
v.e.v., determining the values of baryon masses [1]-[11]. (Here q = u,d are the fields of 
u, d-quarks). In this case O c v _ = 1 and in accord with (11), (13) we have 

n = -(0 | qq | 0) = (1.65 ± 0.15) x 10~ 2 GeV 3 (at 1 GeV) [10] (14) 

(The integration over SU(2) subgroup position in SU(3) colour group as well as anti- 
instanton contribution are included in the definition of n.) The anomous dimension of 
quark condensate is equal to 4/9. According to (14) n has the same anomalous dimension. 
The size p of the zero-mode wave function can be found by calculation in the framework of 
our model of the v.e.v. 

- \ n 

-9(0 | ^—G^J> | 0) = ml(0 | qq | 0). (15) 
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The parameter m 2 , is equal to [13]: m\ = 0.8 GeV 2 at 1 GeV. The m\ anomalous dimension 
is equal to -14/27. Working in the SU{2) colour group, substitute \ n by r a (a = 1, 2, 3) and 



take for G a the instanton field 



4 p 2 

G%{X, P) = -^ (a , 2+p 2)2 » ( 16 ) 



where the parameter r\ ailv were defined by 't Hooft [14] (see also [5]). The substitution of 
(12) and (16) into (11) gives after simple algebra 



1 1 

2 p 

Therefore, 



n— = rn^n. (17) 



= 0.79 GeV' 1 = 0.256fm(at 1 GeV). (18) 



V2m 

This value of p is close to ones, used in the delute instanton gas or instanton liquid models. 

We are now in a position to calculate less well known quantities - the magnetic suscep- 
tibilities of quark condensate, induced by external constant electromagnetic field. 
The dimension 3 quark condensate magnetic susceptibility is defined by [15]: 

(0 | qa^q | 0) F = e q x(0 | qq | 0)F^ U , q = u,d, (19) 

where quarks are considered as moving in external constant weak electromagnetic field F^ v 
and e q is the charge of quark q in units of proton charge (the proton charge e is included 
in the definition of F^). The left-hand side of (19) violates chirality, so it is convenient to 
separate explicitly the factor (0 | qq | 0) in the right-hand side. It was demonstrated in [15] 
that (0 | qa^ u q \ 0)p is proportional to the charge e q of the quark q. A universal constant x 
is called the quark condensate magnetic susceptibility. 

Let us determine the value of x m our approach. For this goal it is necessary to consider 
Eq.3 in the presence of external constant electromagnetic field F^ and to find the first order 
in F^ u correction to zero mode solution (12). This can be easily done by representing ip as 

il>(x,p)=ij) Q (x,p)+il) 1 (x,p), (20) 

where ipo is given by (12) and ipi represents the proportional to F^ v correction. Substitute 
(20) in Eq.3 added by the term of interaction with electromagnetic field, neglect ipi in this 
term and solve the remaining equation for ipi(x, p)). The result is: 

1 / lx 2 \ 

ipi(x,p) = -^e q r] a ^a a F^x 2 (l + -— jip (x,p), (21) 

where a a are Pauli matrices. The matrix element ip + a^ v ip appears to be equal: 

= ~e q F^+x 2 {l + \^)^o- (22) 

(The properties of r\ ailv symbols [14], [5] were exploited.) The v.e.v. (19) in the Minkowski 
space-time is given by: 

(0 | ?(vV I 0) F = e q F^n± J dV + ^ f p2)3 (23) 
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(The normalization condition (13) for ip (x,p) was used.) It is convenient to express n 
through quark condensate by (14), use the notation x 2 = r 2 , where r is the radius- vector in 
4-dimensional space. Then according to (19) we have: 

The integral (24) is quadratically divergent at large r. So, the cut-off R is introduced. Its 
value can be estimated in following way. The volume occupied by one zero-mode in 3- 
dimensional space is approximately equal to \ jn (the volume of the Wigner-Seitz cell). So, 
for cut-off radius square R 2 in four-dimensions we put 

A / \ 2/3 

R2 = 1(A\ =7 .92GeV- 2 
3 V4W 

where the factor 4/3 corresponds to transition from 3 to 4 dimensions. The calculation of 
the integral (24) at the values of parameters p (18) and R 2 (25) gives 

X = -3.52GeV~ 2 (26) 

The quark condensate magnetic susceptibility was previously calculated by QCD sum rule 
method [16]-[19] and expressed through the masses and coupling constants of mesonic reso- 
nances. The recent results are: 

x(lGeV) = -3.15 ± 0.3GeV~ 2 [18]; *(1 GeV) = -2.85 ± 0.5GeV~ 2 [19] (27) 

(The earlier results, obtained by the same method, were: x(0-5 GeV) =—5.7 GeV~ 2 [16] and 
x(lGeV) = 

-4.4±0.4 GeV -2 [17].) The anomalous dimension of x is equal to -16/27. It was accounted 
in [16]-[19], but not in the presented above calculation. (In some of these papers, the a s - 
corrections and continuum contribution, were also accounted.) One can believe, that the 
value (26) refer to 1 GeV, because the value of quark condensate (14) refer to this scale and 
also because the scale 1 GeV is a typical scale, where, on the one hand, the zero-modes and 
quark condensates are quite important (see, e.g. [12]) and, on the other, the instanton gas 
model is valid [9]. Since the integral is quadratically divergent it is hard to estimate the 
accuracy of (26). I guess, that it is not worse, than 30-50%. In the limit of this error the 
result (26) is in an agreement with those found in phenomenological approaches. 

Turn now to quark condensate magnetic susceptibilities of dimension 5, k and £ defined 
in Ref.[15] 

9(0 I q\^ n G; u q | 0) F = e q KF^(0 | qq | 0), (28) 

-iffW<0 I STs^G^g | 0)f = e ? £iV(0 | qq | 0) (29) 

Perform first the calculation of k. In this case the expression of ipi(x, p) (21) must be 
multiplyed by the additional factor: \r h G b where G b is given by (16) and the indices p, v 
in (21) are changed to A, a. In the further calculation it will be taken into account, that xo 
in (12) corresponds to total spin-colour isospin J = and consequently 

<tt X o = -3xo, v a T b Xo = S ab Xo- (30) 
5 



(24) 



(25) 



In the relation 

VbtivVbXa = S^\S ua — 5fj, a 5 v x + E^uXo- (31) 

the last term drops out after summation of zero-modes from instanton and anti-instanton 
configuration. The final result for k is: 



k — — [ u 2 du- —( 1 + ^-u] = — i 

J (u+ l) 4 V 2 J 2 



. , 1N 13 1 11 
\n(z+l)- — + — — + -- 



6 z + 1 2(2 + l) 2 3(^ + l) E 



, (32) 



where z = R 2 /p 2 = 12.7. Numerically, we have: 

re = -0.26 (33) 

The calculation of £ is very similar to those of k and the result is 

f = 2k = -0.52 (34) 

The values of k and £ only logarithmically depend on the cut-off. But unfortunately the 
logarithm in (32) is not very large and its main part is compensated by the term —13/6, 
appearing in (32). So, the accuracy of (33), (34) can be estimated as about 30%. The 
phenomenological determination of 5-dimensional quark condensate magnetic susceptibilities 
was performed by Kogan and Wyler [20] along the same lines, as it was done in [16], [17]. No 
anomalous dimensions were accounted. The results of [20] are: 

k = -0.34 ±0.1, £ = -0.74 ±0.2 (35) 

As can be seen, they are in a good agreement with (33), (34). The 5-dimensional quark 
condensate magnetic susceptibilities play a remarkable role in determination of A-hyperon 
magnetic moment [21]. 

I conclude this part of my talk. It was argued, that chiral symmetry violation in QCD 
arises due to zero-mode solution of Dirac equation for massless quark in arbitrary gluon field. 
The model is proposed similar to delute instanton gas model, in which the zero-mode solution 
is the same as in the field of instanton. The parameters of the model: the density of zero- 
modes and their size are determined from the values of quark and quark-gluon condensates. 
In the framework of this model the values of quark condensates magnetic susceptibilities 
of dimensions 3 and 5 were calculated in agreement with ones found by QCD sum rules 
method using the properties of hadronic spectrum. The agreement of these two approaches 
gives the strong argument in favour that the instantons are the only source of chirality 
violating condensates in QCD. 

Turn now to the temperature dependence of quark condensate, considered as an order 
parameter. It is expected that quark condensate vanishes at high temperatures and the chiral 
symmetry is restored. Two possibilities are discussed: the second order phase transition 
and the crossover. In our approach the normalization condition (13) takes place at any 
temperature [22]. So, the temperature dependence of quark condensate reduces to 

a(T) = (0\qq\) T = n(T) (36) 

One may expect that n{T) vanishes only in the case, when the quantum number, curried by 
instanton is vanishing. This quantum number is the topological charge, which is temperature 
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independent. So, in our approach n(T) never vanishes and the phase transition is of the 
type of crossover. The order parameter temperature dependence of the second order phase 
transition near the critical point T c is smeared by fluctuations [23] (The same statement 
refers, surely, to the crossover.) The fluctuation are determined by long wave oscillations of 
the fields. In QCD they are given by small frequency gluonic field G™^, which is a constant. 
We have: 

~ (0 I G 2 | 0)T c V , (37) 

where Aa(T c ) is of the order of magnitude of the variation of ot(T) near the crossover critical 
point. It is evident, that the left hand side should vanish, if T c would be zero. I guess, that 
the factor, corresponding to this circumstance is T c 2 . The factor p 6 is added for dimensional 
reasons. 

All the discussion of the a{T) temperature dependence is based on implicit assumption, 
that this object has a physical meaning at finite T. It is unclear, however, how it can be 
defined. I did not succeed to formulate the gedanken-experiments, in which this object could 
be measured at finite T. Therefore, may be, such an object has no physical sence at T ^ 0. 
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Program grant RUP2-2961-MO-09. I acknowledge the support of the European Community- 
Research Infrastructure Integrating Activity "Study of Strongly Interacting Matter" under 
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